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THE BCS GAP EQUATION FOR SPIN-POLARIZED 

FERMIONS 

ABRAHAM FREIJI, CHRISTIAN HAINZL, AND ROBERT SEIRINGER 

^^ , Abstract. We study the BCS gap equation for a Fermi gas with un- 

Cn ' equal population of spin-up and spin-down states. For cosh(5^/T) < 2, 

with T the temperature and 5^ the chemical potential difference, the 
question of existence of non-trivial solutions can be reduced to spectral 
properties of a linear operator, similar to the unpolarized case studied 
^^ , previously in [7l[8l|9]. For cosh(5^/r) > 2 the phase diagram is more 

complicated, however. We derive upper and lower bounds for the critical 
(~| ' temperature, and study their behavior in the small coupling limit. 

1. Introduction 

Spin-polarized fermionic systems have attracted a lot of interest in recent 
years, and have been intensely studied both experimentally and theoretically. 
^ ' They are of relevance in various areas of physics, ranging from cold gases to 

^^ . nuclear stars. 

^ I The goal of our present paper is to study the BCS gap equation for such 

f^ ■ imbalanced systems. Experiments |171ll8j have demonstrated a significantly 

richer phase diagram than in the balanced, unpolarized case. The possibility 

f^ ' of pairing in polarized systems was first pointed out in [16], and since then 

has been intensively studied in the literature, mostly under the assumption 
of a contact interaction among the fermions. We refer to p] and [3] for 
recent reviews on this subject. 

^ , Our work is a continuation of recent papers [3 El El [TO] where the BCS 

^ '_ gap equation in the balanced case for systems with general pair interaction 

potential V was investigated. In particular, a criterion for potentials giving 
rise to superfluidity was derived. The form of the interaction in actual 
physical systems can be quite general, and hence it is important to keep 
V as general as possible. In the following we recall these results before we 
present the main results in the imbalanced case. 



C^ 



1.1. Balanced fermionic systems. Consider a gas consisting of neutral 
spin ^ fermions. The kinetic energy of these particles is described by the non- 
relativistic Schrodinger operator, and their interaction by a pair potential 
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which we write for convenience as 2\V, with A being a couphng parameter. 
According to Bardeen, Cooper and SchriefFer (BCS) [2] a superfluid state of 
the system is characterized by the existence of a non-trivial solution of the 
gap equation 



at some temperature T > 0, with E{p) = \/(p^ — /u)^ + |A(p)p. Here, fi gM. 
is the chemical potential and V{p) = (27r)~'^'^ f^^ V{x)e''^P^dx denotes the 
Fourier transform of V. The function A(p) is the order parameter which is 
closely related to the wavefunction of the Cooper pairs. 

The detailed investigation of Eq. I^A^ in[71|8l[9lll0l[IIl[l2l[l3]is based 
on the observation in [8] that the existence of a non trivial solution of (jl.ip 
is equivalent to the existence of a negative eigenvalue of the linear operator 
Kri-iV) + XV{x) on ^2(1^3)^ ^^^^^g 

Kt{p) - 



xann — Typ — 



The critical temperature Tc{\V) is thus defined by the equation 

inf spec {Kt, + A^) = . (1.2) 

Since Kt{p) is pointwise monotone in T, this defines Tc uniquely. For T < 
Tc{\V) the equation (jl.ip has a non-trivial solution, while for T > Tc{XV) 
it has none. 

For /u > 0, Kt{p) attains its minimum on the sphere 0^ = {p '■ P^ = m}. 
In the small coupling limit A — )• 0, the study of Tc reduces to the study of 
an effective operator V^ : L'^{Q^) — )• L'^{Q^), given by 

(M(^)=(^^-^j^ V{p-q)u{q)du;{q), (1.3) 

with duj denoting the Lebesgue measure on the sphere. The lowest eigenvalue 
e^ < of V^ determines the leading order of Tc{XV) as A — )■ |7j. In fact, 
in [3 [9] it is shown that 

Tc = ^l f^^— + o{l)] e"/(2^^^) as A ^ , (1.4) 



2^ 



fJ-J 



where 7 ^ 0.577 denotes Euler's constant, and 6^ = 7rAe^/(2y^) -|- 0(A 
where the 0{X'^) term is determined by another effective operator oi L'^{Q 
The parameter b^ can be interpreted as an effective scattering length (in 
second order Born approximation), which reduces to the usual scattering 
length of the interaction potential 2XV in the low density limit, i.e., as 
//^O [lOj. 



the bcs gap equation for spin-polarized fermions 3 

2. Main results 

We consider a gas of spin ^ particles at temperature T > 0, interacting 
via a pair interaction potential 2XV. We assume the two spin states are 
unequally populated, but have equal masses. We introduce two separate 
chemical potentials, n+ and ;U_, for the spin up and spin down particles, 
respectively. We denote the average chemical potential by /2 = (^_|_ + ^_)/2, 
and half their difference by 

(>ii = 2 ' ^ ' 

which we assume to be non-negative without loss of generality. 

Let 7+(p) and 7-(p) denote the momentum distributions for the spin-up 
and spin-down fermions respectively. In addition, let a denote the Cooper 
pair wave-function. The BCS functional J-t is defined as 

^T{l+n-,a) = - {p'^ - fi+)-f+{p)dp + - {p^ - fi-)-f-{p)dp 

+ xf |a(j;)|V(x)dx-^S(7+,7_,a), (2.2) 



(2.3) 



where the entropy S 


is 




5(7+, 7_, a) 


= - [ Trc4[T{p)lnT{p)]dp, 

JR3 


with 

T{p) = 


fl+ip) 




\a{p) 


a{p) \ 
l-{p) -a{p) 
-a{p) l-7+(p) 

l-7-{p)/ 



(2.4) 



The functions 7-1- and a are constrained by demanding that < r(p) < 1 as 
an operator on C^, for all p E M'^. 

The functional J-'t can be heuristically derived [151 [8] by starting from a 
genuine many-body Hamiltonian and making three steps of simplifications. 
First, one considers only quasi- free (or generalized Hartree-Fock [l]) states. 
Second, one ignores the exchange and direct terms in the interaction energy. 
Finally, one assumes translation invariance and works out the energy per 
unit volume. The specific ansatz of the states in the form of F follows from 
assuming, in addition, that the total spin points in the direction of the 
applied field (quantified by 5^). 

We will show below that on an appropriate domain the functional J^t 
attains its minimum and the corresponding Euler-Lagrange equation, ex- 
pressed in terms of the order parameter A = —2XV * a, takes the form 



where E(p) = \J {jp' — /i)^ + |A(p)p. For convenience we define the convolu- 
tion with a factor (27r)^^/^ in front, i.e., f*g{p) = (27r)^^/^ J^-j, f{p—q)g{g)dq. 
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In terms of 



2T J ^ ""^^^^ \^ 2T 

the gap equation can conveniently be expressed as A = — AV^ * [A/K^^ ). 
Define also 

x>0 ^ 

The basis of our analysis is the following theorem. 

THEOREM 1. Let V G L^/^^M^) he real-valued, fl eR, T > and 6,, > 0. 

(i) If the linear operator Kj, g {—iV) + XV{x) has a negative eigenvalue, 

then the functional Tt attains its minimum at some a ^ 0. 
(ii) If J-T attains its minimum at some a 7^ 0, then A{p) = —2XV * a{p) 

satisfied the BCS gap equation ()2.5p . 
(iii) If the BCS equation has a non-trivial solution, then the linear oper- 
ator KT^s^i—i^) + \V{x) has a negative eigenvalue. 

The proof of Theorem [1] is given in Section [3l 

Remark 1. We shall show in Corollary [T] that for all pairs {6^,T) with 
cosh((5^/T) < 2, the function K^^ is pointwise monotone in A and hence 

This means that statements (i) — (iii) in Theorem [1] are actually equivalent, 
exactly as in the balanced case [8l Thm. 1]. Hence we have obtained a 
linear characterization for the existence of a non-trivial solution of the non- 
linear BCS gap equation. This does not hold if cosh{6^/T) > 2, however, as 
discussed below. 

2.1. Critical temperatures. Since we do not have equivalence of the state- 
ments (i) — (iii) in Theorem[T]in case cosh(5^/r) > 2, we do not have a simple 
definition of the critical temperature as in (jl.2p in this case. With the aid 
of the operators Kj.^ and Kj^^s^, we can define curves T^ and T^ in the 
{5f,,T) plane, given by 

TliXV) : {{5,„T) I inf spec(i^f ,^ + XV) = 0}, (2.7) 

T^iXV) : {{6^,T) I inf spec(i?T,v + XV) = 0}. (2.8) 

Below the curve Tj one has inf spec (ivT^ ^ + XV) < 0, and hence The- 
orem [T]^i) implies that the functional J-t has a global minimum at some 
a 7^ 0. Above the curve T^ , however. Theorem [1]^ iii) implies that the gap 
equation does not have a non-trivial solution. Consequently the minimum 
of J-'t is attained for a = 0, the normal state. See Fig. [TJ 

In the region between Tj and T^ , there can exist non-trivial solutions to 
the gap equation even if a vanishes identically for the minimizer of J-'t- In 
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fact it was shown in [5j by numerical calculations, using a contact interaction 
potential, that there is a parameter regime where the gap equation has a 
solution but the corresponding energy is higher than that of the normal 
state. 

Using the methods developed in [3 [9] we are able to evaluate T^ (AF) and 
T^ {XV) in the small coupling limit A — t- 0, and obtain analogous formulas 
as in the balanced case in (II. 4p . With Tc denoting the critical temperature 
in the balanced case (at /i = /i), we shall show that 



■^# 



Tr,e 



-^*(S./Tf) 



(2.9) 



as A — )• 0, with ^ standing for either i or o, and the k* are two explicit non- 
negative functions (defined in (j4.9p and (j4.29p . respectively), with k#(0) = 
0, of course. We refer to Theorems [2] and [3] in Section S] for the precise 
statements. The resulting curves are plotted in Figured! 



5„/T=cosh"'(2) 



VT=1.91 




Figure 1 . The curves T^ , Tf and T^ (in increasing order) , 
in units of T^ the critical temperature in the balanced case. 
In region I the system is in a superfluid phase. In region IV, 
there are no non-trivial solutions to the BCS gap equation, 
and the system is the normal phase. In region III we prove, 
under additional assumptions on V, that the system is in a 
normal phase, even though non-trivial solutions to the gap 
equation might exist. 



2.1.1. A more detailed phase diagram. Under additional assumptions on the 
interaction potential V we are able to derive a sharper upper bound on the 
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critical temperature, at least for small coupling. Let T^ be the curve defined 

by 

THXV) = {i5^,T) : inf inf spec(i^^, + XV) = 0}. (2.10) 

Note the difference with Tg in ()2.8p . where the minimum is taken pointwise 
for every fixed p in KJ^ ^ (p) . Recall also from Remark [1] that KJ^ g is 
not monotone in y if cosh {5^/T) > 2, so the minimum of the spectrum of 
^TS "*" -^^ ^^ ^'-'^ necessarily attained at y = 0. 

In terms of a the BCS gap equation (j2.5p can conveniently be written as 

^i^^,^+Ay)a = 0. (2.11) 

Let us assume that Kf^g -\-XV has as its lowest eigenvalue. For minimizers 

1 /J. 

of the functional J-'t, this is the case, for instance, of V is negative, since in 
this case a is negative and hence the ground state, by a Perron-Frobenius 
type argument. 

In this case, we will show that A can be treated as a constant in the small 
coupling limit, with its value given by the one it attains on the Fermi surface. 
As a consequence, we shall show that if {6^,T) lies above Tf , the system 
is in a normal phase (in the sense that a = for the minimizer of J-t), 
even though the BCS gap equation can have non-trivial solutions. Hence 
the sharper upper bound T^ on the critical temperature for superfluidity 
holds in place of T^ ; see Fig. [TJ The precise formulation of our results is 
given in Theorem [4] in Section [H 

Numerical calculations in [5] (see, in particular. Fig. 2 there) predict that 
between the curves Tj and T^ the gap equation has at least two solutions. 
Moreover, there exists another dividing line inside this region separating 
the parameter regimes where the BCS functional is minimized by a = or 
a 7^ 0, i.e., the normal and the superfluid regime. 

2.1.2. A toy model in one dimension. For the purpose of illustration let us 
now consider the BCS functional for a one-dimensional system, where the 
particles interact via a contact potential of the form V{x) = —g5{x) with 
5 > 0. In this case the gap equation takes a very simple form, and the order 
parameter A is constant. In fact, the BCS gap equation in this case is 

1 1 /■ 1 , , . 

-dp. (2.12) 



9 ^TT^K^^g^ip) 

The critical temperatures Tg and T^ are given implicitly via the equations 

A (2.13) 

dp \ . (2.14) 




- ' / 


' d^ 


2vr 7iR 
1 


f ' 


A 27r ^ 


i^k^Ap) 



THE BCS GAP EQUATION FOR SPIN-POLARIZED FERMIONS 7 

The curve Tg encloses the region where the gap equation ()2.12p has ex- 
actly one solution. Above T^ there are no solutions, while in-between (j2.12p 
has exactly two solutions. Numerically this is easy to see, in fact. If one 
plots 2^ /ig T^A^ I \ dp as a function of A, one observes that the graph crosses 

1/g either once, twice, or not at all. 

3. Preliminaries and proof of Theorem [T] 

We start by specifying the precise domain of definition of the BCS func- 
tions J't in dil. 

DEFINITION 1. Let V denote the set of functions (7+,7_,a), with 
7+,7_ G L1(M3^(1 +p^)dp), < 7+(p),7_(p) < 1, and a £ H^(R^,dx), 
satisfying |a(p)|^ < 7+(p)(l -7-(p)) and |a(p)P < 7-(p)(l -7+(p))- 

For (7+,7_,a) G V, the corresponding T{p) in (j2.4p satisfies < F < 1 
as an operator on C^ for all p € M^. Moreover, all the terms in Tt are 
well-defined under our assumptions on V. 

Recall that fl = {/j.^ + fJ-~)/2 denotes the average chemical potential and 
6^ = (/x+ — H-)/2. It is easy to see that in the absence of a potential 
V, Tt is minimized on D by the normal state given by a = 0, 7*J_(p) = 
[eW-f^y^^^ + l]-i and 7° (p) = [eW'fi)+^^^ + l]-i. 

We start with the observation that J^t attains a minimum on D. 

PROPOSITION 1. There exists a minimizer of Tt in D. 

The proof is analogous to [8, Proposition 2], and we skip it for simplicity. 

Let us now come to the proof of our main Theorem [TJ We start with 
some preliminaries. Diagonalizing F allows to rewrite the entropy (|2.3p in 
the form 

5(F) = — / [{r + w) ln(r + w) + {r — w) ln(r — w) 

+ {s + w) ln(s + w) + {s — w) ln(s — w)]dp, (3.1) 
where the functions s, r and w are given by 

r= 2(1 + ^^+-^-)' ^= 2^^"^+ + ^-)' '^ = -\/(l-7+ -7-)^+4|ap. 

Lemma 1. Let (74-,7_,a) £ D be a minimizer of Tt for < T < 00. 
Then, for a.e. p G M^, 

(V * a)ip) = -^a{p) ifriw) + fs{w)) (3.2) 



2 2 \s^ -w 



^. = ^^^^ = ^1- :T^)' (3-3) 



(p2 -fl) = ^(l-^_^- l-){fr{w) + fs{w)), (3.4) 



where fa{b) = lln^,. 
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Proof. LetA,:={pGR^:{\aip)\^--/+il--f-)){\aip)\^-^-il--f+))>e} 
for some e > 0. Let g S //^(M^)nL^(M'^) be such that the Fourier transform 
^ of gf is supported in A,:. Then for small enough \t\, (7_|.,7_, a + tg) S V. 
By the Lebesgue dominated convergence theorem, it follows that 



7 /" 

— J"T(7+,7->a + *5)t=o =2 K / a{x)g{x)V{x)da 



+ |k I a{p)g{p)ifr{w) + fs{w))dp. (3.5) 

Here ^ denotes the real part. Also, for g real, (7+ + tg, 7_, a) € P for small 
|i|, and we get 



—J^Th+ + tg,-f-,a)t=o = - 
^- , (7++7--I) 



and, similarly, 
d 



dt 



^t{i+,1- +tg,a) 



t=o 



T 

+ 2 



(7+ + 7- - 1) 



{p^ - fj-+)g{p)dp 



(/rH + /sH)+ In 



(/ - f^~)9ip)dp 



{fr{w)+fs{w))-ln 



w 



w 



w 



w 



dp (3.6) 



dp . (3.7) 



Since (74.,7_,a) minimize Ft by assumption, the expressions in (J3.5p . (j3.6 
and (j3.7p vanish. This implies that (j3.2p holds, and also that 



(/-/x_) = - 



^^ ^+ ^-\fr{w) + f,{w)) + \n(-^ 



w 



2\ 1 



w^ 



(3.8) 



(/ - /X+) = - 



(1-74 



^^{friw) + fsiw))-ln'''^ 



w 



2\ n 



W 



(3.9) 



for a.e. p € ^e- As in [8j we can argue that the measure of i? := M \Ue>o^e is 
zero. Subtracting and adding Eqs. (j3.8p and ()3.9p implies (j3.3p and (|3.4p . D 



Proof of Theorem [Jl We consider the case where T > in the proof, and 
leave the analogous case T = to the reader. We proceed similarly to ^ 
Theorem 1]. To show (i), note that 



KP 



T,5, 



ip) 



2(p2 _ fl) 



P - fJ' 



tanh 



2T 



+ tanh 



p'^-fi-5^, 
2T 



1 



7+ -7° 



(3.10) 



where 7^ and 7^ are the momentum distributions in the normal state defined 
in the beginning of this section. We shall show that if the normal state a = 
ininimizes J-t, then inf spec {K^ ^ + Ay) > 0. 
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Since it is true that for every g G C(f , ||(1 - 7+ - 7_) + (ifi') ||oo < 1 
for small enough |t|, it follows from the Lebesgue dominated convergence 
theorem that -^r^J^Til+^l-^tg) exists for small t, and equals 



^^Trijljttg) 



t=o 



2A / \gix)\^V{x)dx + T / \g{p)\ 



ln(fej)+ln(fej; 



dp 



at t = 0. Here r' 







^(1+7^ 



7- 



2(1 — 7+ — 7-) • It is easy to see that 



In 



7.0 _ ^0 



+ ln 






(1-7^-7°) 

^(1 - 7^^ +7°) and w^ 

2(p2 _ fl) 



2A 



I \g{x)\^V{x)dx + 2 1 kIs,{pMp)\^^P- 



and hence 

^J-T(7+,7-,t5)t=o 

If the normal state minimizes Tt then clearly Sr^J'Til+il- 1^9)1=0 ^ 
(keeping in mind that it is also stationary). This implies that {g, {K^ ^ + 
^y)g) ^ for all g £ C^ , and hence proves the claim. 

Next we shall show (ii). For a minimizer (7+,7_,a) of J-y, we can define 
A via 



A{p) 



p"^ - H 



1 - 7+ - 7- 



-a{p) . 



(3.11) 



Setting further E{p) = \/(p^ — /i)^ + |A(p)P and recalling that 

2w= x/(l-7+-7_)2 + 4|d|2 
we obtain that 

E{p) = 2w{p) 

If we further use ()3.4p we get 



P -^ 



l-7+(p)-7-(p) 



E{p)=W-{fr{w) + fs{w)). 

It then follows from the definition of A(p) and E{p) that 

A(p) 



Noting that 



E{p) 

ifriw) + fs{w)) 



w{p) = a{p). 

(r + w)(s + w) 
In ■ 



(3.12) 
(3.13) 



oj {r — w){s — w)^ 
we obtain from (j3.3p and the definition of E(p) that 



E{p) + S, 



In 



r + It; 



w 



E{p) - 5, 



In 



s + w 



w 



10 A. FREIJI, C. HAINZL, AND R. SEIRINGER 

Since it is also true that 

1 - 1 + —— - 1 + ^— = 2u;, 3.14 

V s — w \ r — w 



we get from (j3.2p and (|3.13p that 



A(p) 
and hence 



\V*a = =^^, (3.15) 



A / 1 

A = -2\V *a = -\V * — 1 



^ I l^r+w -^^ 



T+W 

s— ui ~ ' r—w 



From the above expressions, we observe that 

E{,p)+s^ r -\- w E{p)~s^ s + w 

e T = ^ e T = . 

s — w r — w 

Therefore we arrive at the BCS gap equation 

A / , /, . ;^kp)+Sji.\~'^ ( ^ g(p)-v 

'e 



A = -Ay*— ( 1- ( 1 + e T \ - ( 1 + e t ) ). (3.16) 



-1^ 



To see (iii), note that the gap equation can be written as 

with a(p) = ^[p)/Kj, p [p). Consequently, 

(a, {kT,6, + XV)a) = (a, (Kt,,5, - ^^^Ja) < (3.17) 

by the definition of -fCr.iS • The inequality is, in fact, strict if A is not 
identically zero, since a and A have the same support. This completes the 
proof of Theorem [TJ D 

COROLLARY 1. Assume that cosh(5^/T) < 2. Then the existence of a 
non-trivial solution of the BCS gap equation is equivalent to K^ g + Ay 
having a negative eigenvalue. 

Proof. The corollary follows directly from Theorem [1] and the fact that 
K^g (p) is pointwise monotone in A for all S^/T E [0, cosh" (2)]. Hence 

K^g = Kt,5^ in this case. To see this one checks that the function 
X !->• (tanh(x+c)+tanh(x— c))/a; is monotone decreasing on M+ if cosh(2c) < 
2. D 
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4. Bounds on the critical temperature 

As explained in Subsect. 12. H the two linear operators in Theorem [1] de- 
termine upper and lower bounds T^ and T^ on the critical temperature. 
In this section we evaluate Tg (XV) and Tg {XV) in the limit of small A. We 
will proceed similarly to [71 [9l [TT] . Note that it follows from Corollary [T] 
that for cosh{5^/T) < 2 the two curves coincide, i.e., T] (XV) = T^ (XV). 
(Compare with Fig. [TJ) 

We start with Tg (A^), defined in (j2.7p . The Birman-Schwinger principle 

implies that a pair (5^,T) € Tg {XV) is characterized by the fact that the 
compact operator 

Bt,s, = X\V\-2j^V-2 (4.1) 

has —1 as its lowest eigenvalue. Here we use the notation 

V{x)^/^ = {sgnV{x))\V{x)\'^/^. 

That —1 has to be the lowest eigenvalue of -Bt,(5^ follows from the monotonic- 
ity of Kj,g in 6^, similarly to [71 Lemma 3.1]. In fact, since the spectrum of 
Bx^Sf^ goes to for large 5^ and is continuous in 6^, an eigenvalue of Bt^s^ 
smaller than —1 would correspond to an eigenvalue —1 for larger 6^, which 
is in contradiction to the monotonicity of Kj, g + XV in 5^ and being the 
lowest eigenvalue for the critical parameter. 

For fl > 0, the operator Bt,s^ becomes singular as {6^,T) — )• (0,0), and 
the key observation is that its singular part is represented by an operator 
V^ : L^Qf,) ^ L\nf,), given by 

where duj is the uniform Lebesgue measure on 0^. We note that the operator 
V^ has already appeared in the literature O [H] . If we assume that V G 
-L^(M^), then V(p) is a bounded continuous function, and hence Vp is a 
Hilbert-Schmidt operator. In fact, V^ is trace class, and its trace equals 



27r^ 



/jg3 V{x)dx. Let 



infspecV^ (4-3) 



denote the infimum of the spectrum of Vn. Since Vn is compact, we have 



e/i<0. 



Let J- : L (R ) — )• L (Qfi) denote the bounded operator which maps 
ip S L^(R^) to the Fourier transform of ip, restricted to the Fermi sphere 
Qp. For V S L^(R^), multiplication by l^p^ is a bounded operator from 
L^(R^) to L^(R^), and therefore J"|y|^/^ is a bounded operator from L'^' 
to L\nf,). 
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Finally, as in [9j, we define Wp, via its quadratic form, to be an operator 
on L^(J7^) such that 



HW^|n)= / d|p|(— '^'' 



dn{mp)\^-mvnp/\p\)\' 



+j^ I dn\^{Vfip/\p\)\'] ■ (4.4) 



Here, {\p\, Q) £ M4. x S^ denote spherical coordinates for p G M^ and (p{p) = 
(27r)~^/^ /^ V{p-q)u{q)duj{q). It was shown in [7] that (|4.4p is well defined, 
despite the apparent singularity of the integral, and that W^ is of Hilbert- 
Schmidt class. 
For A > 0, let 

Bp = X-^V-^-X^^Wi,, (4.5) 



'^-''2^^'^ 2/2' 



and let ^(A) denote its ground state energy, 

^(A) = inf spec^S^ . (4-6) 

We note that if e^ < 0, then ^(A) < for small A. If the (normalized) 
eigenfunction u G L'^{Qp) corresponding to e^ is unique, then 

Q{X) = {u\Bfi\u)+o{X^). (4.7) 

In the degenerate case, this formula holds if one chooses u to be the eigen- 
function of Vfj, that yields the largest value (u|VV^|u) among all such (nor- 
malized) eigenfunctions. 

4.1. Evaluation of T^ . In this subsection we will derive the asymptotic 
behavior (j2.9p of Tj {XV) for small A, which we present in the following 
theorem. 

THEOREM 2. Let V G L^{R^) D L^/'^{R^) be real-valued and ft > 0. 
Assume that e^ defined in {-j-^ ^^ strictly negative, and let q{X) he defined 
as in ( (^.g| ). Then any sequence of pairs ((^^(A),T(A)) on the curve Tg {XV) 
satisfies 

1'-?. ('" f + viw - '''<*'"'^0 = ' - ^ - '"'*/" <'^** 

where 7 ~ 0.5772 is Euler's constant and 
K\t) = -/ ^ — + 7/ 7 In-. (4.9 

Note that n'{t) = ln(i) -h 7 - ln(7r/2) -h o(l) for large t, while k;^(0) = 0. 
We can rewrite Eq. (|4.8p in the form 
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If Tc denotes the critical temperature at 5^ = 0, defined in Subsect. ll.lt this 
can also be written as 



K'^^J=ln^^J+o(l), (4.10) 

which agrees with ()2.9p and is plotted in Figure [TJ The region denoted by I 
defines the region where the system is guaranteed to be in a superfluid phase. 
The asymptotic behavior of k} implies that the curve hits the vertical axis 
at 1 and the horizontal axis at ne"^ /2 ~ 0.88. 



Proof. Let 



m{T,6^) = ^[ I ^ l^]dp, (4.11) 

and let M(T, 6^) be the operator 

M(r, 5^) = j^ - m{T, 5^,)FF* , (4.12) 

with J^ defined in the paragraph after Eq. (14. Sh . Following \T, Lemma 2] it is 
not difficult to see that V^/^M{T,6^)\V\^/^ is bounded in Hilbert-Schmidt 
norm uniformly in (T, 6^). 

Suppose ip is an eigenstate of the operator K^^ + XV corresponding 

to the lowest eigenvalue 0. By the Birman-Schwinger principle, 4> = V^ip 
satisfies 

with i?T,5 defined in (|4.1|) . Moreover, —1 is the lowest eigenvalue of i?T,5 ; 
as argued above. 

Since V^''^M{T,6^)\V\^''^ is bounded uniformly in {T,6^), the operator 
1 + \V^''^M{T,5^)\V\^''^ is invertible for A small enough. Therefore 

1 + My\ ^ T^^^ = 1 + A|y|^ (m(T, 5,,)FF* + M{T, 5,,)) v'^ (4.13) 



'T,5^ 






= Dh+ Xm{T, 6^)D-'V2T*T\V 

where D = 1 + XV^M{T,6^)\V\^ . The fact that -1 is an eigenvalue of (fiTJl 
yields 

D ll + Xm{T,6u) ^ -V^T*T\V\^] = 0. (4.14) 

V 1 + XV2M{T,6^)\V\2 J 

Since D is invertible, this is equivalent to the operator 

Xm{T,6.)T\V\^/^ -—4 ---n:,V^/^T* (4.15) 

having an eigenvalue — 1 (using that AB and BA are isospectral if they are 
compact). Note that the latter operator is acting on L^(f]p). 
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By expanding (1 + XV'2 M{T,6fj_)\V\^)~^ in a Neumann series we obtain 
the following implicit characterization of the corresponding temperature: 

Xm{T, 5^)mi spec (J^ {V - XVM{T, 5^)V + ©(A^)) T*) = -1. (4.16) 

Note that J^VJ'* = y/p.Vp,, which was defined in (|4.2p . It follows that to 
leading order 



1 

A-s>o ^ ' ^^ inf spec-Fyj-"* \/^^, 



hm^ Xm{T, 6^) = --^— — ^.^^^ = -^^-. (4.17) 



To derive the second order correction we employ W^, which was defined in 
(j4.4p . It follows from (j4.16p . first order perturbation theory, the fact that 
J-VJ-* is compact and that inf spec J-"yj-"* < by assumption, that 

"'^^' ^^^ " X{u\TVT*\u) - X^u\TVM{T,6f,)VT*\u) + 0{X^) ' ^^'^^^^ 

where u is the normalized eigenfunction corresponding to the lowest eigen- 
value of TVJ-*. If u is degenerate, we choose the u that minimizes the A^ 
term in the denominator of (I4.18P among all such eigenfunctions. Eq. (j4.18p 
represents an implicit equation for T| . Since J-VM{T,6^)VJ-* is uni- 
formly bounded and T — )• as A — )• 0, we have to evaluate the limit of 
{u\TVM{T,5^)VT*\u) as {5^,T) -^ 0. For this purpose, let c^ = VT*u. 
Then 

{u\TVM{T,6^)VT*\u) 

— g — —\0{p)fdp-m{T,6f,) / \ip{p)\'^ duj{p) 

( K^n^ [l^b)l' - \^{Vnpl\p\)?] + ^\^{V^P/\P\)\ 




Recall the definition of Kj.^ (p) in (|2.6p . For a fixed Sfj_/T, it converges to 

\p'^ — p.\ as T — 7- 0. Using the fact that the spherical average of |c^(p)p is 
Lipschitz continuous (see [9]), we can interchange the limit and the radial 
integral over \p\, and hence obtain 

lim{u\TVM{T,6f,)VT*\u) = {u\Wp\u) . (4.20) 

We have thus shown that for (5^,T) E Tg 

hm (m{T, 6^) + — —^ -^-—) = . (4.21) 

A^o V mf spec (A VA* Vp. - A^ W^) / 

It remains to compute m{T,5^). 
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Lemma 2. As {T,6f,) -^ (0,0), 

= -^(l.i^ + 7-2 + li.^-Kl(Vr) + o(l)). (4.22) 

Proof. We split the integral into two parts according to whether p^ < fi or 
P^ > /2. By changing variables from p^ — p, to —t and t, respectively, we see 
that m(T, 5^) equals 



J_ f'' // _ 1 1 \ V^^^ _ _J__\ 

+ ^r((l '- '- )^^^^ —]dt 

(4.23) 
To simplify the notation, let us introduce the function 

^o(t) = f 1 It-tt It-yt) ■ (4-24) 

We can rewrite ()4.23p as 

™(T.«^i/;(T„,o:^-^)* 






dt 



1 r^ / 1 1 \ 1 r^ To(t) 

' ' + ,_ dt + ^= / ^^^dt . (4.25) 



2fJ, Jq \^/^ + t y/n-tj \^ Jo t 

Using the dominated convergence theorem, it is easy to see that 



oo 



Moreover, 



hm / {To{t)^^^^^-^L=]dt 



I r^ f I 1 , , 

+ ,_ . \dt 



Again by the dominated convergence theorem, the second integral becomes 
J-r-/^ + ^^^-^^dt=^(ln4-ln(l + ^) + V2-2 
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in the limit T, 5^ -> 0. With c = 5^/T, using the fact that [1 + e"]''^ + [1 + 
e"^]""*^ = 1, we can rewrite the last term in ()4.25p as 



1 


\r^u ' M ix 1 /* V ' ' ]dx] 


v^ 


Jo x^l + e" l + e^e^)"^'^ ' Jq x \l + e-" 1 + e-^e^ ^ ""^ 


Note that 




[^U ^ Mdx - ' [^^-'~^dx 


Jq X yi + e" 1 + e^^e^ J 1 + e^ 7o 2; 




1 fr I — e~^ dx 


l + e" Jq 1 + e^e^ x 


Using integration by parts, 


/ ^ dx = \u^(l e-~^'^) / ln(x)e-^dx. 
Jo X T \ J Jq 


Moreover, since 




- / e-''ln(x)dx = 7 (4.2 



Jo 
(Euler's constant), we conclude that 



lim I / dx — In — = 7 • 



T,5^^0 \ Jq X T 



The same argument applies with c replaced by — c, and hence 
1 rTo(t) , 1 / , /x 8 :, , 

Combining all the terms gives the statement of the lemma. D 

Theorem [2] follows from (fOT]) . (jM}, and (g^!]). D 

4.2. Evaluation of T? . For x > and c > 0, consider the function 

/(a;,c) 



tanh ^ + tanh ^ 



For each c there is a 6(c) such that x 1— )• /(x,c) attains its minimum at 6(c) 
and / is monotone increasing for all x > 6(c). Recall from Corollary [1] that 
6(c) = for c < cosh^^(2) . We have 

K^. =Tf (V(p2-/2)2 + |A|Vr, 6,/T 



Therefore, K^g is monotone in A whenever '^ j,^' > b {5^/T), and we have 
that 

_ ' Kl,^ for^>6(VT) 

(4.27) 



" ^ '" ' 2r/(6(vr),5^/r) for^<6(Vr). 
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THEOREM 3. Let V G L1(M3) n L3/2(m3) he real-valued and fi > 0. 
Assume that ep defined in 14. 3\ ) is strictly negative, and let g{X) be defined 
as in i4.6^ . Then any sequence of pairs ((^^(A),T(A)) on the curve Tg {XV) 
satisfies 

a ('" f + 27m> - «°(«./^)) = 2 - 7 - 1.. I . (4,28) 

where 7 ~ 0.5772 is Euler's constant and 

1 f"" 1 - p-^ dr 1 /""^ 1 - p"^ dr 

K^ic) = ^— / \ — + —— / — (4.29) 

b 



+ ( 1 - e-' 



In 6 — / ln(x)e ^ dx — - — - — In — , 

Jo ^ ' V{b,c) 2' 

with b = 6(c) defined above. 

Recall that 6(c) = for c < cosh" (2), and hence k°(c) = k'(c) in this 
case. For larger c they differ, however, as Figure [1] shows. 

Proof. The proof works analogously to Tg . The only difference is that 
m(r, 5^) now has to be replaced by 

fh{T,S,) = ^f ( \ ^^]dp, (4.30) 

and the operator M(T, 6^) by 

M(r, (^^) = -^ \ - m(r, 5^)-F^* . (4.31) 

'f^r,5^(p) 

The result is that 

lim ( m(T, 6^) H 7 ^ r | = (4.32) 

^^o\ ^ ' ""' infspec(A^V^-A2W^)y ^ ^ 

for {5^,T) £ T^ . What remains is to calculate rh{T,5p). 
Lemma 3. As {T,6^) -^ (0,0), 

m(r, 5^) = ^ An ^ + 7 - 2 + In ^ - K°(5^/r) + o{l)] . (4.33) 

Proof. For simplicity we will abbreviate 6 {5fj_/T) by 6 throughout the proof. 
First note that from ()4.27p it follows that m{T, 6^) equals 



1 



A-KfJ, 



^->AK.sAp) pA'^''^hj^<Xf{h,5,lT) p^H^ 



For given 6 > 0, 



(4.34) 



ri™ ATTfi l\2^<, \2Tf{b, 6^/T) p^ ) '^^ 2^^/(6, 5^/T) ' 
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and the result obviously also holds for 6 = 0, where both sides are zero. 

As in the proof of Lemma O we split the first integral in (14.34p into two 
parts according to p^ < fi and p^ > ji, and change variables from p^ — fl 
to —t and t, respectively. Introducing again the function Tq in (j4.24p . the 
integral equals 






dt 



1 r( 1 ^^\„,^^rMi,,. (4.35) 



2a^ JTb \Vf^ + t V/" - V y/i^ JTb t 

The first integral converges to jl~^''^ In (l + \/2) as T, 5^ — >• 0. The second 
integral gives -y^ (— ln(l + \/2) + -v/2 + In 4 — 2) and the third integral is 

\jll]i in this limit. 

To evaluate the last term in (j4.35p . we proceed as in Lemma [2l with the 
obvious modifications. The result is that 

Collecting all the terms proves the lemma. D 

This concludes the proof of Theorem O D 



4.3. Evaluation of Tf and a more detailed phase diagram. Our final 
goal is derive a sharper upper bound on the critical temperature which sepa- 
rates the phases where the minimizer of the BCS functional has a vanishing 
or a non- vanishing a. We are able to do this under additional assumptions 
on the interaction potential V ^ which, in particular, imply that the BCS 
minimizer a is the ground state of the operator K^ ^ + )\V . 

Recah the definition of Tf in (f2T0]l . 

THEOREM 4. Let V he a radial function in L^/'^{R^)nL'^{R^), with V <0 
and V{0) < 0. Assume that e^ = infspecV^ < 0, and let q{X) be defined 
as in ( f^.Tp . Then any sequence of pairs ((^^(A),T(A)) on the curve Tf {XV) 
satisfies 
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where K^{t) = mfrf>o C(iirf)) with 

({t,d)=- f e-"" In {x/d)dx+ f e'"" In (l + ^Jl - {d/xfjdx (4.37) 
1 /•°° l-e"^ dx 



i. + e-' Jd i + e-^"- Vx^ - d^ 

Invr . 



1 + e-* 7rf 1 + e^-* Vx2 - d2 

Moreover, for any sequence {6^{X),T{X)) that lies strictly above the curve 
Tf , in the sense that 



t/ie -BC^ 5ap equation (j2.5p cannot have a solution for small enough A wii/i 
i/ie property that the operator K^^ + AT/ /las as lowest eigenvalue. In 
particular, the normal state minimizes Ft in this case. 

We note that linid^o ^(t, 0) = K}{t). Moreover, K^{t) = K}{t) for t < 
cosh" (2) ~ 1.32. Numericany, one can check that this equaUty holds even 
on the larger interval t S [0, 1.91]. (Compare with Fig. [TJ) 

The last statement of the theorem concerning the fact that the normal 
state minimizes J-t follows from a Perron- Frobenius argument, as explained 
in Subsect. 12.1.11 It implies that a non- vanishing and minimizing a is nec- 
essarily the ground state of K^^ + W , with eigenvalue zero. 

Proof. By the previous arguments we obtain for {5^.,T) G Tf the following 
behavior in the small coupling limit: 

1 



where 



\\Tn\mi5a,T) + —- ,^ ^^, — — — = , (4.39) 






4vr/i y Jrs I -fCf . ^ (p) p2 



For the first part of the theorem it remains to evaluate m (5^,T), which is 
done in the following lemma. 

Lemma 4. In the limit (T, (5^) — )• (0,0), 

m(<5^,T) = -^(ln^ + 7-2 + ln^-Kg(VT)+o(l)) . (4.40) 

Proof. As in the proof of Lemma [2l we can rewrite rn{5^,T) as 

m((^^,r)=maxl((5^,r,y) (4.41) 

y 
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where 



/(*„, T, y) ^± 1°° (t. (t) ^j^ - ^) * (4.42) 

1 /"A X it) 

1 fP f 1 ^ 1 \ ^^ ^ 1 fP Ty{t) 



, , ^ , ^ ,dt + —= / , ^' ' dt 

2/^70 VV^ + * VfJ'-iJ VJ^ Jo y/t^ + 2/2 



and 



T,(t) 



1+e T 1+e T 



It is clear that the maximum is attained for a y that goes to zero as 
{T,5^) — )• (0,0). Using the Lebesgue dominated convergence theorem, 
one observes that in the hmit T,y,6^ — t- the first integral converges to 
/i-V2 1^ (1 + ^) ^ the second becomes /i^^/^ ^_ ^^(^ _ i) + ^ + in 4 - 2) , 

and the third \/2/fi. 

It remains to compute the last integral in (j4.42p . Changing variables to 
X = yt^ + y^/T, we can rewrite it as 



1 



1±^=^±^1^ dx = I + II := (4.43) 

d Vx"^ - d"^ 



y/x'^ — d^ Jd \/x'^ — d^ 



where c = ^ and d = U^. Now / in (j4.43p equals 



1 + e^ id 



1 



Vx2 - (i2 Vx2-d2 Vl + e^ 



c+c 



dx , (4.44) 



and similarly for //, replacing c by — c. For the second term in the integrand, 
we can simply replace the upper integration boundary by oo as T — )• 0. To 
evaluate the integral of the first term, we integrate by parts and obtain 



^<PHn/T?\ 1, / /i , J^2 , iniT\^ 



1 _ g-v^^+W^r ) In ( ^ + ^^2 + {fx/TY ) - ( 1 - e-'^) Ind (4.45) 



ry/dUiP/Tf _ 

/ e "^ In ( X + \/x2 — d'^j dx . 
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The latter integral converges as T — )• 0. Proceeding in the same way with 
//, we have thus shown that 



^J■ 



V^T7 



dt (4.46) 



— /"CxD 

ln('^) +ln2- {l-e-'^\\nd- / e"^ In (x + \/x2 - dA dx + o(l) 
1 /""^ 1-e-^ dx 1 /""^ l-e"^ dx 



l + e^ Jd 1 + e^'+^ Va;2 - d'^ 1 + e"^ 7^ 1 + e^"^ Vx2 - d"^ 
Combining all the terms, keeping in mind (j4.26p . we arrive at (|4.40p . D 

We now turn to the proof of the second part of Theorem HI The order 
parameter A(p) is assumed to satisfy the BCS gap equation. Under the 
assumption that a is the ground state of K^^ + W we can use the Birman- 
Schwinger principle to show, as in [9', Lemma 4], that 

A{p) = -g{X)l f V{p-q)du:{q)+XMp)], (4-47) 

with ||z^a||oo ^ C uniformly in A, and g{X) a normalization constant de- 
termined by the gap equation. Eq. (j4.47p can be derived using standard 
perturbation theory applied to the operator 

Xm(T,d.,AmV\^/^ ^ — -^ r--T77^V^/'^J^* (4.48) 

where 



m{T,5f„A) = —^ { -TT^—r^ - — ] dp , 
^TTf^ Jr3 \K^^s^{p) p^J 

and M is defined as in ()4.12p . with K^ g replaced by K^g and m(r, (5^) 
by m{T,dfj,,A). The eigenvector corresponding to the lowest eigenvalue —1 
of this operator is, to leading order, given by the lowest eigenvector of Vp,, 
which is the constant function u{p) = J-^ . In fact, the radial symmetry 
of V implies that the eigenstates of Vp, are the spherical harmonics, and the 
fact that V < forces the ground state of Vp to have a fixed sign, hence 
to be constant. The lowest eigenvector of (|4.48p is, therefore, (/> = u + X^\, 
with ^x uniformly bounded in L'^{Q.p) norm. With the aid of the Birman- 
Schwinger principle, a can be recovered from (/> via a = cV^i'^F*(l) for some 
normalization constant c. Hence we have 

A{p) = 2K^^s,{pMp) = -2XV^{p) = -2Xc\V[^(^{p) , 

which implies (j4.47p . 

Using the Lipschitz continuity of the main term in (j4.47p . we can argue 
as in [Hi Thm. 2] that for small A the function m {6^, T, A) is determined by 
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the value of A(p) at the Fermi surface p^ = jl, i.e., 

m i5f„ T,A)=m {6^,T, A(V^)) + o(l) . 
This imphes further that 

hm (m(5„,r,A(V^)) + — 7T-^, r^JT77T)=0. (4.49) 

A^oV ^' ^ infspec(AV^V^-A2W;i)y ^ ^ 

Since, by definition, 

m{5^,T,A{^/fl)) <m{5^,T), 

it follows from (j4.39p that (|4.49p cannot be satisfied for (5^,T) outside of 
Tf , i.e., for such {6^,T) such that 

lim ( m {S^, T) + — — -i- TsTTTT ) < . (4.50) 

A^o V ^ mfspec(A^V;2 - A^W^i)/ 

This completes the proof of the theorem. D 
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